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Starting from a phase-field description of the isothermal solidification of a dilute binary alloy, 
we establish a model where capillary waves of the solidification front interact with the diffusive 
concentration field of the solute. The model does not rely on the sharp-interface assumption, and 
includes non-equilibrium effects, relevant in the rapid- growth regime. In many applications it can 
be evaluated analytically, culminating in the appearance of an instability which, interfering with 
the Mullins-Sekerka instability, is similar to that, found by Cahn in grain-boundary motion. 



PACS numbers: 

INTRODUCTION 

Crystal growth from an under-cooled fluid phase is 
frequently described by a diffusion equation for heat or 
particle density, complemented by boundary conditions 
at a moving sharp interface between the solid and fluid 
phases. As explained in the review [1] by Langer, one of 
the boundary conditions comprises the conservation of 
the density field in terms of the related flux and source 
terms at the interface. A second boundary condition is 
the Gibbs-Thomson relation which is adequate in cases 
of local thermodynamic equilibrium. For applications in 
the rapid-growth regime Aziz and Boettinger [2] have 
extended the Gibbs-Thomson relation by including non- 
equilibrium effects which they derived in an intermediate 
step from an atomistic picture of the initially extended 
interface region. 

One of the most exciting achievements of the sharp- 
interface approach is the explanation of dendritic growth 
in the diffusion-limited regime, opening up at low under- 
cooling. The initiation of this process is described by the 
Mullins-Sekerka instability [3] whereas the cumbersome 
route to an analytic calculation of the fully-developed 
three-dimensional dendrite by Brener [J is reviewed in 
Ref. [5J. Another exciting case of pattern formation in 
the kinetics-limited rapid-growth regime is the periodic 
formation of layers with alternating homogeneous and 
dendritic micro-structures in dilute binary alloys [6j. A 
crucial point in explanations of this effect by Carrard 
et al. [5J, and by Karma and Sarkissian [7] is the non- 
monotonous dependence of the interface temperature on 
the growth velocity, observed in Ref. [2]. 

A fundamentally different access to a theory of crystal 
growth rests on the use of a phase-field order parameter 



which allows a continuous transition between the solid 
and fluid phases. This approach is closely related to that, 
used by Halperin, Hohenberg, and Ma [H] for studying 
dynamic critical behavior in their model C , and, in the 
context of crystal-growth processes, has been described 
by Collins and Levine [9], by Caginalp and Fife [10], and 
by Langer pj] . Phase-field models for binary alloys have 
been established by Wheeler et al. [T2], and later by Kim 
et al. [13 . In the article [TT] by Langer it is emphasized 
that one reason for promoting the phase-field approach is 
the hope, that similar to the identification of universality 
classes for dynamic critical phenomena, one may obtain 
a fundamental understanding of various pattern-forming 
mechanisms in solidification problems. A quite different 
reason for the interest in the phase-field approach is that 
it presents a convenient basis for simulations, avoiding 
the costly interface-tracking procedure, necessary in the 
sharp-interface description. 

Within a phase-field approach Lowen et al. p3] have 
investigated the long-time crossover behavior between 
the diffusion- and the kinetics-limited regime. One of 
their results is the appearance of meta-stable branches 
where the steady-state motion of a planar solidification 
front can occur in the diffusion-limited regime. These 
branches are part of the trajectories which show the non- 
monotonous dependence of the interface temperature on 
velocity, later found in Ref. [2]. In a phase-field model 
of rapid solidification Ahmad et al. [H] have considered 
the effects of solute trapping and solute drag, and found 
partial agreements with the results, obtained in Ref. 0. 
A good example for the computational use of the phase- 
field approach is the work of Karma and Rapell [16] on 
a quantitative description of dendritic growth. In their 
treatment the above-mentioned non-monotonic velocity 
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dependence of the interface temperature is noticed by 
the appearance of a negative kinetic coefficient, joining 
the term, linear in velocity. Since such an anomaly does 
not arise in the sharp-interface description of dendritic 
growth, the authors invented a compensation device for 
the unwanted term. 

Among the numerous papers on derivations of sharp- 
interface descriptions from a phase-field model the work 
by Elder et al. [U] is probably the most elaborate one. 
Their approach avoids the assumption of a zero interface 
width, but instead uses the products of this width with 
the interface curvature, and with the growth rate divided 
by the diffusion constant, as small expansion parameters. 
Since, however, the second parameter is of order one in 
the rapid-growth regime, applications to this regime are 
excluded. The result for the Gibbs-Thomson relation in 
the low-velocity regime contains a kinetic-under-cooling 
term where, as in Ref . [16] , the kinetic coefficient can, in 
some parameter range, become negative, an effect which 
regrettably has not been scrutinized by the authors. 

In the present paper we advertise a model description 
where, retaining a finite interface width, the interface 
position is used as basic field variable, in addition to a 
bulk-diffusion field. Although the energy density is a 
feasible example of a bulk field, we will predominantly 
consider the isothermal solidification of a dilute binary 
alloy, in which this role is taken by the concentration of 
the solute component. An advantage of our approach is 
that it remains valid in the rapid-growth regime which 
can be seen in a thorough derivation from a phase-field 
model. Our description can, however, also be derived in 
a self-contained way from first principles, in the first line 
from the observation that the presence of an interface 
breaks Euclidean symmetry, and, therefore, implies the 
existence of Goldstone modes [TH]- These are capillary 
waves of the solidification front which interact with the 
bulk-diffusion mode. The interaction kernel as well as the 
diffusion coefficient can be freely chosen in our approach 
which, therefore, can be applied to a sizable variety of 
model systems. 

For a subset of models with a location-independent 
diffusion constant, we have derived a universal form of 
the dispersion relation of interface eigenmodes. These 
modes determine all possible morphological instabilities 
of the solidification front, and in all cases we encounter 
the Mullins-Sekerka instability. It is the only instability 
which survives the sharp-interface limit in our scheme. 
As soon as we allow a finite interface width, we become 
aware of an additional instability which features a finite 
amplification rate already for a planar perturbation, a 
behavior, previously discovered by Cahn [T!5] in grain- 
boundary motion. The nature of this instability is closely 
related to the non-monotonic velocity dependence of the 
interface temperature. In our approach the shift of this 
temperature from the melting point of the solvent enters 
as a driving force which we also evaluated analytically 



for several models. Although these subjects are, like the 
distantly-related Corriel-Sekerka instability [20], mostly 
discussed in connection with the rapid-growth regime, 
they also affect the low-velocity behavior, as illustrated 
by the previously mentioned sign problem of the kinetic 
coefficient. 

One of the toy models, investigated in the frame of 
our capillary- wave approach, allows a surprisingly simple 
analysis of some attributes of the solidification process. 
In this model both input-functions, the solute-interface 
interaction and the solute-diffusion coefficient, are taken 
to linearly interpolate between the solid and fluid bulk 
phases, assuming a finite interface width. The effect of 
solute trapping is measured by the partition coefficient, 
for which we found an expression, coinciding with that, 
derived by Aziz and Kaplan |21) . up to some re-scaling 
of the involved characteristic velocity. We, furthermore, 
have established a simple analytic relation between the 
driving force and the growth rate which reflects the non- 
monotonous behavior of the interface temperature. The 
minima of all driving-force trajectories are connected by 
a kinetic spinodal line which we also have determined 
analytically, finding agreement with an analogous line in 
the non-isothermal growth of a one-component crystal, 
discussed by Umantsev [35]. The image of this spinodal 
line in the temperature-concentration phase diagram is 
located between the solidus line and a static spinodal 
line, established by Baker and Cahn [23]. Within our 
model this provides an answer to an issue concerning this 
matter, discussed by Hillcrt in Ref. [24 . 

In our concluding discussion we present an expression 
for the entropy production in the steady-state growth 
of a planar solidification front. The result demonstrates 
that the appearance of a negative kinetic coefficient, also 
in our approach, is not in conflict with basic principles 
of linear irreversible thermodynamics. Also included in 
our discussion are some estimates for the relevant model 
parameters which determine the range of validity of our 
approach. We, finally, emphasize the flexibility of our 
description, concerning generalizations and the inclusion 
of additional field variables. 



PHASE-FIELD DESCRIPTION 

In order to derive our capillary-wave model from a 
standard phase-field model, we initially describe the 
isothermal solidification of a dilute binary alloy by a 
phase field $(r, t) for the solvent, and a concentration 
field C(r,t) of the solute component. In terms of these 
field variables the effective Hamiltonian of our model 
reads 
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FIG. 1: Double-parabola potential W($) with equilibrium 
phase-field values $ = +1 in the liquid, and <£> = — 1 in the 
solid phase. 



FIG. 2: Temperature-concentration phase diagram, showing 
the liquidus and solidus lines Tl{C),Ts{C) which enclose the 
two-phase region. The values Cl and Cs define the miscibility 
gap at the temperature Tl(Cl) = Tg(Cs). Also shown is the 
equilibrium spinodal line Tb(C) according to Baker and Cahn. 



H 



d A rl - 



rt2 



y(V$) 2 + W($) 



k(T) 



[C -U(<f>,T)]' 



(1) 



Here, W($) is a double-well potential with degenerate 
minima at $ = +1 and $ = — 1, representing the liquid 
and solid equilibrium phases of the solvent, fn the present 
paper we will exclusively use the explicit form 



W($) = 6(-$) - ($ + l) 2 + 6($) - ($ - l) 2 , (2) 

shown in Fig. 1. Contrary to that, the potential W(<£>, T) 
is adaptable, up to the relations 



U{-1,T) = C S (T) , U(+1,T) = C L {T) 



(3) 



where Cs (T) and Cl (T) are the solute concentrations in 
the solid and liquid phases at temperature T, depicted in 
the temperature-concentration phase diagram in Fig. 2. 




C 



FIG. 3: Double-tangent construction, applied to the free- 
energy densities fs{C,T) and fh{C,T). The intersection 
point Co(T) defines the Tb-line in Fig. 2. 



by the double-tangent construction, visualized in Fig. 3. 
In this procedure the terms, involving the equilibrium 
chemical potential /j,e(T), cancel. The intersection point 
in Fig. 3 at 



The values ^ follow from the mean-field free-energy 
densities 



fs(C,T) 
h(C,T) 



2 

k(T) 



[C - U(-1,T)} 2 + fi E (T)C , 

[C -C/(+l, T)} 2 + fi E (T)C , (4) 



Co(T) 



C S (T) + C L (T) 



(5) 



borders the regions where the two bulk phases can exist 
in a meta-stable state. In Fig. 2 the inverse function 
T (C) defines a kind of spinodal line, promoted by Baker 
and Cahn [33]. Within the often- used approximation of 
a constant miscibility gap this line has the form 
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T (C)= T ^ + ™. (6) 




1- 






The equilibrium chemical potential [i E (T) in Eqs. (|4|) 
obeys the Clausius-Clapeyron equation 










T dflE L (7) 
M dT ~ AC{T) K ' 
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where Tm is the melting temperature of the pure solvent, 
L the latent heat, and 




-1- 







AC(T) = C L (T)-Cs(T) , 



(8) 



is the miscibility gap. Following Langer [T], we can use 
Eq. ([7]) to extract an expression for k(T) by forming the 
total derivative of the chemical potential 



FIG. 4: Single-kink phase-field profile in thermal equilibrium, 
following from the potential W( ( I > ). 



H L (C, T) = k(T)[C - C L (T)} + p E (T) 



(9) 



with respect to temperature, taken at C — Cl{T). The 
result reads 



«(T) 



fdC L 



L 



\dT j T M AC(T) 



(10) 



where, as already in Eq. ([T]), only leading terms in AC 
have been taken into account. We later will realize that 
the assumption of a small miscibility gap is a crucial 
assumption in the derivation of an interface description 
from a phase-field model. 
The equilibrium conditions 



5J±_ SH 
5$ ' SC 

have the single-kink solution 



(11) 



r — 



5H 
5$ 



d t C = V-X>($)V 



T 

1 SH 

K SC 



(14) 



where T is the phase-field relaxation rate, and T>($) the 
diffusion coefficient of the solute atoms. The relations 



V(+1)=D L , V(-l) 



D, 



(15) 



allow in general different values Dl and D$ in the liquid 
and solid phases. Furthermore, 



J r ($) = F[9(-$)($ + 1) - 9($)($ - 1)] 



(16) 



is a driving force which is operating near the interface, 
and, for F > 0, supports the solid at the expense of 
the fluid phase, thus activating the solidification process. 



The form ( 16 ) can be viewed as arising from a potential 



$ E (z) = e(-2)[-l+exp(V0] 
+e(z)[l-exp (-*/£)] 



(12) 



describing the solid-liquid phase-field profile, displayed in 
Fig. 4, and an attached solute-concentration profile 



C e {z)=U[$e(z)] 



(13) 



which leads to identify 2£ with the width of the interface, 
and the parameter a in the Hamiltonian ([!]) with the 
surface tension. 

For the dynamics of the system we adopt the equations 
of motion 



w F ($) = e(-$) 
+ e($) 



1 + —F I ($ + l) 2 - —F 



a 



1 — —F I ($ — 1] 
a 



(17) 



F 



replacing W($) in the Hamiltonian |l]). A convenience 
of this modelling is that it remains meaningful even close 
to the apparent spinodal point at F = noted in 

Rcf. 25J. We finally mention that Langevin forces, 
representing thermal-noise effects, have been ignored in 
Eqs. (14), because we are primarily interested in the 



non-equilibrium effects, induced by the force T . 

A dimensionless form of our model equations can be 
established by the mappings 
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AC 



Di 

— > r 

(C - C S ) - 



■t->t 



(18) 



U 



F 



where the shift of the field C has been made with regard 
to the case of an under-cooling at constant concentration 



Cs- The equations of motion (14) then assume the form 



P 



d t C = V -VV{C-U) . 
Here, two independent parameters 



F + 1 {C-U) 



dU 



(19) 
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P : 



Vc 



(20) 



appear, the latter expressed in terms of a crystallization 
and a diffusion velocity, 



v c = re , v D 



Dl 



In the following we are going to reduce Eqs. (19) to 



an interface description, valid up to growth rates of the 
order of Vjy, and keeping a finite interface width 2£. This 
will be done, adopting the specific forms ([2| and ( 16 ) for 
W(<I>) and but successively using different choices 

for the parameters 7 and p, and for the functions U{&) 
and £>($). 



SINGLE-COMPONENT SOLIDIFICATION 

We first consider the case 7 = which offers to describe 
the solidification of the pure solvent by the remainder 

1 



- d t <5> = V 2 $ - 6(-$)(l + F)($ + 1) 

p 

-e($)(i-F)($-i) 



(22) 



of Eqs. ( 19 ). In this scheme the propagation of a planar 



solidification front with constant velocity 

V = vV D (23) 
in z-direction, is represented by the single-kink solution 



(24) 



in the co-moving frame where Eq. ( 22 ) reduces to 



- q(-0(<s>f + 1) - e(0(*f - 1) = 



(25) 



P 



& F + F[Q(-Q($ F + 1) - 0(C)(*F - 1) 



Assuming that $f is an odd function of £, both sides of 



Eq. ( 25 ) vanish separately. As a consequence one finds 



MO = e(-c)[-i + exp(c)] 

+e(C)[l-exp(-C)] 



(26) 



which is identical to the equilibrium phase-field profile 
$e(0- Moreover, Eq. (25) fixes the relation between 



the growth rate v and the driving force F in the form 



v 
P 



(27) 



The simple results (26) and (27) are due to the specific 
choices ^ and( 16 ) of our model. Especially, the former 
one will allow us to largely copy the derivation of an 
interface description near thermal equilibrium, presented 
inRef. [26]. 

In order to study small perturbations of the steady- 



(21) state solution (26), we linearize Eq. (22) in 



(p(x,z- vt,t) = $(r,f) - $ F (z-vt) . (28) 
This leads to the equation 



dt-d* + Sl)<p(x,C,t) = 



(29) 



involving the two-dimensional Laplacian 

d 2 = V 2 - dl , (30) 
and the one-dimensional Schrodinger-like operator 



n = -a 2 



Fd c - 26(0 + 1 + F[e(-C) - 6(C)] • (31) 



By taking the derivative of Eq. (25) with respect to £ 



one finds that $^(C) is an eigenfunction of f2 with zero 
eigenvalue. It is identical to the bound state of the delta 



potential in Eq. (31 1, and also follows from translational 



symmetry of the system in ^-direction. The excited state, 
closest to this ground state, can be found by application 
of the pseudo-gauge transformation 



' X P ! +7T C I exp [ -— c 



-dl 2 «y(o + (i + f ) 2 ®(-0 + (1 - f ) ' 



(32) 



9(0 
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The operator (32 1 has a band of eigenstates 



MO = e(-c)exp(o 
+e(C) 



(33) 



K — 2 

cos (kQ H — sin (fcO 

k 



above the ground state which are parametrized by the 
wave number k. The related eigenvalues read 



e{k) 



(34) 



and are identical to those of the operator fl. The second 



identity in Eqs. ( 34 ) implies that the band of eigenvalues 



e(k) is separated by a gap Ae > 1/4 from the ground- 
state eigenvalue e = 0, provided F < 1. This constraint 
excludes the regime beyond the ghostly spinodal point 



of the potential (17 1, and, in view of the result (27), is 



equivalent to the statement V < Ve- 



in Eq. ( 29 ) solutions of the form 



^(x, C> *) = $ 'f(C) exp (i q • x - et) 



(35) 



represent a band of soft modes which have been identified 
in Ref. [TS] as Goldstone modes due to broken Euclidean 
symmetry in space. They are well separated from the 



hard modes, arising from the excited states (331, and 
physically are over-damped capillary waves. 

In order to implement the interface position R as a 
collective coordinate in the co-moving frame, we follow 
Ref. [57] where any point Q = (x, £) near the interface 
is represented in the form 



Q = R(s,i) + uN(s,t) . 



(36) 



Here, N is a unit vector, normal to the interface, as shown 
in Fig. 5, and 



u = u(x, £, t) , s = {s , s } = s(x, C, t) (37) 

define mappings to a normal coordinate, and to a set of 
curvilinear coordinates within the interface, respectively. 
With the notation di = d/ds l the rate of any quantity 



*(8 > ti J t) = *(x,C,t) (38) 
at constant x and £ can be written in the form 



fltf = (3 u tf ) dtu + (S^) d t P + dt* 



(39) 



N(s,t) 




FIG. 5: Definition of the local curvilinear coordinates {s, u} 
where Q = {x, £} is some position close to the point R(s, t) 
at the interface, and N(s, t) is the normal vector at R(s,t). 



The rates d t u and d t s l are given by the relations 

d t u = -N • d t R , (40) 
( 9ij - u K i3 ) d t y = - {8 t R + u d t N) ■ djR 



which follow from Eq. ( 36 1 , and involve the metric and 



extrinsic curvature tensors of the interface, 
gij = (9iR) • djR , 



(41) 



As shown in Ref. [27], the metric tensor of the three- 



dimensional embedding space, induced by Eq. (36), is 
block diagonal with G uu = 1, and with the 2-dimensional 
sub-matrix 



Gi 



(g ik - uK ik ) g (gi 2 - uK t j ) 



(42) 



where gijg^ k = Sf. This form is needed to evaluate the 
Laplace-Beltrami representation 



diVGGP'dj + d u VGd u ) * , 



again with G tJ G^ k = 6$, and G = Det{G i3 }. 
We next decompose ^(r,^) in the form 

$(r,i) = <f> F (u)+r](s,u,t) , 
along with the Fadeev-Popov condition [55] 



(43) 



+00 



du& F (u) r]{s,u,i) = 



(44) 



(45) 
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This condition ensures that rj(s, u, t) is a pure hard-mode 
field, and implicitly defines R(s, t) as a collective field 
variable. Since excitation of the contribution 77 requires 
to overcome the gap A e in the eigenvalue spectrum of 



the operator (31 1, 77 may be treated as a perturbation in 
Eqs. (22 1. An expansion of this equation to linear order 



in 77 yields 



isothermal conditions. We mention that the derivation of 
this result neither relies on a sharp-interface description 
nor on a restriction to the low-velocity regime. 

For practical calculations it is convenient, to evaluate 
Eq. (48) in a Monge representation with respect to the 
planar steady-state motion. In this representation the 
interface-position and normal vectors are given by 



tf F (u) i-dtu 



V 2 u 



P 



0(77) =0 



(46) 



where the zeroth-order term derives from Eqs. (39) and 
(43), applied to *B(u) = $>p(u). All other terms of zcroth 



order cancel due to the identity (25), with ( replaced by 
u, and v by the normal velocity component vjy = vN z . 
If, by multiplication with & F (u) and integration over u, 
Eq. (46) is projected onto the soft-mode subspace, the 



terms linear in r\ either drop out due to Eq. ( 45 1 , or they 



carry pre- factors of the type ( 40 1 or ( 41 ) , and, due to 
this, are of higher order in our basically hydrodynamic 
approach. 

In view of Eqs. (40) and (43) this procedure leads to 
the result 



- N • <9 t R = 

P 



du [& F {u)} 2 Tr 



K 



-uK 



(47) 



adopting matrix notations for the metric and curvature 
tensors. The singularity, emerging in the integrand of 
this expression, if u approaches the smallest curvature 
radius of the interface, is cured by the factor [$^(u)] 2 
which decays on the scale of the interface thickness. As 
discussed in Ref. [27], this allows to expand the singular 
part in Eq. (47) in powers of u, giving rise to moments 
of the weight [$' F (u)] 2 . Following this line, one finds, up 
to higher-order curvature corrections, 



P 



N • <9 t R = Tr [K • g" 1 ] = K 



(48) 



where K is the local mean curvature of the interface. 



We now add the result ( 48 1 as a perturbation to the 



R(s,t) = {s, h(s,t)} , 

N(b,4) = —±={-dh(a,t),l} 

Here, h(s,t) is a local height variable, d = 
in-plane nabla operator, and 



g(s,t) = l + [dh(s, t)f 



(51) 

(d 1 ,d 2 ) the 
(52) 



is the Jacobian of the interface metric. The metric and 
extrinsic-curvature tensors read in Monge representation 



9ij 



(dih){djh) , K. K 



V9 



d l d J h . (53) 



Implementation of these expressions into Eq. ( 48 ) leads 
to the differential equation 



— d t h = n5 VJ d l —d 1 h . 
V9 Vfl 3 



(54) 



For a stability analysis of the uniform reference motion 



(25) it is sufficient to linearize this equation in h(s,t). 



The result is a two-dimensional diffusion equation with 
a diffusion constant p which, within this model, implies 
stability of the planar interface morphology. 



In terms of physical units, the nonlinear equation (54) 
can be written in the equivalent form 



— d t h = a SHd 



(55) 



relation ([27}, and, for an under-cooling from the melting involving an effective Hamiltonian 



temperature Tm to some temperature T < Tm, assume 



F = L 



T 



T 



(49) 



and an Onsager coefficient 



(56) 



After going back to physical units via Eqs. (18), (20), 



and (21), we then arrive at the Gibbs-Thomson relation 



A 



re 

a 



(57) 



T M -T 



T, 



= -aK 



M 



(V + N • d t R) 



(50) 



including a kinetic under-cooling term. It locally relates 
the growth rate to the curvature of the interface under 



This representation can be read as a dynamical version of 
the drumhead model which, complemented by Langevin 
forces, has been established and analyzed with regard to 
critical fluctuations in Ref. [25] . 
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BINARY-ALLOY SOLIDIFICATION 

The inclusion of a finite amount of solute particles is 
most easily accomplished, if in the model equations ([!]) 
and (14 1 we choose 



u($) = Cl + Cs + Cl ~ Cs $ 

K 1 2 2 

in accordance with the constraints ([3]), and 



(58) 



Elimination of C F from Eqs. (62) and (64 1 leads to 



the third-order differential equation 



(d c +v) 



-^-6(-C)(l + F + 7 )($ F + l) 



-6(0(1 -F + 7 )($ F -$i) 



which implies the matching conditions 



--y&F (66) 



= D L = D S = D 



(59) 



Adopting this together with the choices ([2| and |l6| ) for 
W and F, the scaled equations of motion ( [19] ) assume 
the explicit form 

= p[V 2 $-6(-$)(l + F)($ + l) 

- 6($)(1 - F)($ - 1) + 7 (C - 1 - $)] , 
d t C = V 2 (C-$) (60) 



where, according to Eqs. (18), C(r,t) means the scaled 
excess concentration with respect to Cs- 
In the one-dimensional stationary case, 



C(r,t) = C F (z-vt)=C F (Q , 



(61) 



the last equation in Eqs. (60) can be integrated once. 



Since, by definition C F {— oo) = 0, one obtains 



C' F 



vC F = <J?'f 



which, due to the behavior $' F (+oo) = C' F (+oo) 
implies the steady-state boundary condition 



C F (+oo) = C F (-oo) = 



(62) 
= 0, 

(63) 



In the equation 



& F -e(-()(l + F + 7 )(<S> F + l) 



-6(0(1 - F + 7 )($ F - $ L ) = -7^ 



(64) 



resulting from the first of the Eqs. (60), the condition 



( 63 1 enforces the identification 



1 - F - 7 

$ F (+oo) = = $ L 

1 — -F + 7 



(65) 



The deviation of this expression from the static value 
$£(+00) = 1 is negligible, if 7 <C 1 — F, as also pointed 
out in Ref. Q3]. 



$ F (-0) = $ F (+0) = , 
$' F (-0) = $' F (+0) , 
$ F (-0) = $ F (+0) + 2 , 



(67) 



including in the first line the definition $ F (0) = of the 



steady-state kink position. Eq. ( 66 ) has the single- kink 
solution 



$ F (0 = 6(-0[-l + ^-exp( a _0] 
+ 6(0[$L + Socxp(/3 

+ B+exp(p+0] (68) 

where a_ , /?o , /3+ are the roots of the cubic equations 



(a + v) 
(P + v) 



a' 



~a- (1 + F) 
P 



P 2 + -f3 - (1 - F) 
P 



(69) 



with a_ > 0,3?(/3 ) < 0, < 0. The amplitudes in 

the solution ( 68 1 , and the connection between _F and v 



are determined by the relations 



A_-1 = B + B+ + *l = 0, 

- £ A) - S+/3+ = , 
A_a 2 _ - B /? 2 - B + f3 2 + = 2 , 



(70) 



following from Eqs. (67) and (68). 

In the limit 7 — > the controlling roots of Eqs. (69) 



approach, due to Eq. (27), the values 



Pa 



(71) 



Accordingly, the term cx Bq in the solution ( 68 1 features 



an exponential long-distance behavior which elucidates 
the non-monotonic phase-field profile, found in Ref. |30j . 
More importantly, however, this behavior endangers the 



suppression effect of the singularity in Eq. (47). A way 



out of the problem derives from the observation that the 
amplitude Bq of the dangerous term obeys the relation 
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10 3 F 



Bo 



Pl + -P a -{l-F) 



P 



(72) 



v _ B v_ 



— 2 l ~ F 
v 1 - F + 7 



which follows from Eqs. (69 1 and (70 1. This relationship 



suggests, and, in fact, forces us to evaluate the model by 
a perturbation expansion in 7. 

In the following we choose to consider the quantity 
v/p = V/Vc as another small quantity. This constitutes 
a noticeable restriction to the growth rates V of materials 
where Vc obeys an Arrhenius law, like in intcrmctallic 
compounds. In dilute metallic alloys, however, particles 
need not overcome an activation barrier in order to form 
a crystalline structure. This has been pointed out by 
Aziz |31j . and is more explicitly demonstrated by the 
molecular-dynamics simulation of the crystallization of a 
Lennard- Jones liquid in Ref. [32] . Consequently, in such 
materials our approximation scheme allows applications 
to the rapid-growth regime where growth rates can be of 
the order of the diffusion velocity Vd- In view of these 
features we now will analyze the solidification process of 
a binary alloy to lowest order of a double-expansion in 7 
and v/p. 

Starting point of the expansion are the relations 



1 + F + j 1 - F + 7 



2 

^ r+00 



$2 



(73) 



P 

Cf(C) 



J — 00 

C dC'$' F (C')exp[ V (C'-C)] 



which directly follow from Eqs. ( 64 1 and ( 62 ) . Neglecting 
terms of order 7 2 and jv/p in the first of these equations 
allows the replacement $f(C) — ► ^s(C)- This leads to 
the expressions 



F(v)-G F (0) = -v-G F (v) , 
P 

C F (C)= f dC'$' B (C')ex P KC'-C)] 

J —OO 

where, in the first equation, the quantity 

+00 



(74) 



g f = -7 / dC^(C)C^(0 

J —00 

/• + 00 

= - jv d([C F (C)] 2 (75) 



A 


/ Y 


V i / 


2 



FIG. 6: The external driving force F, considered as a function 
of the growth rate v according to Eq. (76 1, for p = 100, and 
for 71 = 0.001 and 72 = 0.01, respectively. For the latter 
value F(v) has a minimum, to the left of which the anomaly 
F'(v) < is visible. 




FIG. 7: The steady-state concentration profile for the velocity 
values vi = 0.001, i>2 = 0.01, and ^3 = 0.9, illustrating the 
solute-trapping effect. 



line of Eq. ( |74| ) we have used the identity 27 — — G F (0). 
This term derives from the contribution (k/2) [[/($)] 2 in 
the Hamiltonian (JlJ which, due to the expression ( [58] ) for 
[/($), favors the solid phase, and, accordingly, acts as an 
internal driving force in the steady-state growth process. 
The second line in Eq. (75) follows from Eqs. (62), and 



(63), and implies that the right-hand side of Eq. (74) is 



positive. 



Insertion of the kink solution (26) into Eqs. (74) and 



( 75 ) leads to the explicit results 



defines the solute drag in line with the definition, used 
by Hillert in Ref. [33]. On the left-hand side in the first 
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F(v) = - + 7 
P 



v + 2 
(^Tlj 1 



- 2 



(76) 



c F (c) = e(-0— exp(c) 

v + 1 



-0(0 



1 



exp (-C) 



v — 1 1 " v 2 — 1 
The first of these equations implies 



exp (-uC) 



F'(0) 



1 



1 

Pc 



Pc 



1 



Accordingly, F(v) has a positive slope at the origin for 
p < Pc but shows the anomalous behavior F'(0) < 
for p > pc, pointed out already in Ref. [14]. In the 
literature this anomaly has repeatedly been encountered 
again, however, without any clarification of its physical 
background. Below, we will demonstrate that the effect 
is due to a strong instability of the system, suggesting 
that the inclusion of random forces in the basic equations 



of motion (14) might become important if F'(v) < 0. 
The results (76) are illustrated in Figs. 6 and 7 in the 
unstable regime p > p c - We mention that the existence 
of an anomaly of the type F'(v) < is also in line with 
the velocity-dependent interface temperature, obtained 
in the approach by Aziz and Boettinger [2] . 



The results ( 76 ) also determine the partition coefficient 



K(V) 



C s 



C s + C F (0)AC 



(78) 



which measures the solute-trapping effect, and for V = 
reduces to the equilibrium value 



K E = K(Q) = 



Cs 
C L 



(79) 



By insertion of the value Cf(0) = l/(t> + 1), following 
from Eq. ( 76 ) , one recovers the form 



ip(x,z-vt,t) = $(r,i) - $>f(z - vt) , 
n(-x,z-vt,t) = C(r,t) - C F (z- vt) . (82) 

The resulting equations can be written in the form 



d t — p d 2 
d 2 d t - d 2 



3 



-PI 
—d 2 — vc^ 



(83) 



(77) 

where, replacing Eq. (|31|), now is defined by 



n = -d 2 c - -u c 



dc-2- 



1 



*'f(0) 

f[9(-c) - e(c)] . 



5(0 + 1 



(84) 



By taking first and second derivatives of Eqs. (64) and 



(62 ) with respect to £, one finds that the matrix operator 



on the right-hand side of Eq. ( 83 ) obeys the relation 



81 



-PI 
~~d 2 — vd^ 



C' 







(85) 



which means that translational symmetry in ^-direction 
again provides an eigenstate with eigenvalue zero. 

A convenient way, to explore the appearance of other 
eigenstates, is to expand the elements of the matrix in 
Eq. (85) in the small parameters 7 and v/p. To leading 
order the resulting eigenvalue equation reads 



ptt E 

d 2 — d 2 — i>c\ 



= pe 



4> 



where the operator Qe is given by 



-d 2 Q + 1 - 26(C) 



(86) 



(87) 



The upper component of Eq. ( 86 ) yields the autonomous 
eigenvalue equation 



K(V) 



k e + v/vg 
1 + v/ys 



(80) 



suggested by Aziz in Ref. [3T]. Remarkably, however, we 
moreover find a reference velocity 



V* D = 



Vd 



which shows a dependence on Ke of the type proposed 
in Ref. [34]. 

In order to examine the stability of the steady-state 
solutions C F , we next expand the exact equations of 



motion ( 60 1 to linear order in the perturbations 



n E ^(0 = ei>(0 

whereas the less interesting lower component in principle 



allows to calculate 0(0- Eq. (88) possesses two classes 



of solutions, corresponding to the ground and scattering 
states of the operator ( [87] ) . The ground-state equation is 
a relict of Eq. (85), implying e = and ip(0 — &e(0- 



(81) The scattering states form a band with 



e(k) = l + k 2 



ipk(0 « 



exp (ikO — 



1 



i + |fcCI 



exp(z|fcCI) 



(89) 



(90) 
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Again we have the situation that the soft modes of the 



Eqs. (95) - (97 1 are the main results of the present 



system, deriving from Eqs. (83) and (85), are separated section. Parallel to the procedure, leading to the Gibbs- 



from the hard modes by a gap which, in leading order of 
our expansion, is given by Ae = 1. Since this value will 
only slightly be shifted in our expansion, we will eliminate 
the hard modes of the system by suitably extending the 



procedure, described below Eq. ( 44 ) 



As a first step we complement the decomposition ( 44 ) 
by splitting C into three contributions, 



C(r,t) 

Here, by definition, ■& is related to 



Cf{u) + c(s, u, t) + $(s, u, t) 
C F (u) + c(x,C,*) + tf(x,C,*) 



T](s,U,t) = 77(x,C,t) 



by the equation 



(d t -vd c ~W 2 )d = -V z fi 



?2~ 



(91) 



(92) 



(93) 



which is a copy of the second line in Eq. ( 60 ) , written in 



the co-moving frame. In view of Eq. (62) the remaining 
equation for c reads 

(d t - vd c - V 2 ) c = -{d t u - V 2 u) C' F - (V 2 u)*b (94) 

where we have, in the spirit of our expansion scheme, 
replaced Q' F by When the solution c of this equation 



is introduced into the first line of Eq. (60), projection 



onto the soft- mode component 3>' E (u) leads to the result 



- N • D t R = K - g 
P 

where the last term appears as a perturbation 



+oo 

(s, t) = — 7 / du &e(u) c(u, s, t) 



(95) 



(96) 



of the drag force (75 ). Here, in addition to the neglected 



higher-order terms, leading to Eq. (48), we have omitted 



terms of order 777. The singularities, arising from the 
terms V 2 u = — 7Y[K/(g — uK)] in Eq. (94), are cured 



by the factor 3>' E (u) in Eq. (96 1 . This justifies to replace 
these terms by —K, just as in Eq. (48). As a result, the 



equation ( 94 ) can be written in the form 



(dt - vd ( - V 



2 ^c 



(N • D t R -K)C' F + K& E (97) 



where the differential operator on the left-hand side can 
be represented in terms of curvilinear coordinates via 



Eqs. (39) - (431 



Thomson relation (50), we are going to add the result 



( 95 ) as a perturbation to the force balance in the first 



line of Eq. (74). For an under-cooling from an initial 



point Tl(Cs) at the liquidus line to some temperature 
T < Tl(Cs) at constant concentration we assume 



F = 4 7 



47 



T L (C S )-T 
Tl(Cs)-T s (C s ) 
T S (C S ) - T 
Tl(C s )-T s (Cs) ' 



(98) 



Here, the pre-factor 47 has been chosen such that, after 
transforming to physical units, we obtain 



Ts(C s ) - T 



Ti 



M 



G F (0) -G(s,t) -trK (99) 



DV C 



(V + N • <9 t R) 



which in a natural way generalizes the previous form ( 50 ) 
of the Gibbs- Thomson relation. The new term, involving 
Gf(0) and the total drag force 



G(s,i) = G F +g(s,<) 



(100) 



acts as an additional under-cooling. A similar effect is 
present in the result for the interface temperature in the 
Aziz-Boettinger approach [5] . We mention that Eq. ( [99] ) 
also provides a new answer to the question in Ref. [35] . 
concerning a Gibbs-Thomson equation, valid for non- 
planar interfaces. 

The rather involved evaluation of Eqs. (95) - (97) is 



considerably simplified, if the Monge representation of 
these equations is only used in linear order in h and c. 
This is, fortunately, adequate for a stability analysis of 
the planar morphology of the interface, and leads to the 
differential equations 



d t h = p [d 2 h 



] , (101) 
d 2 ) c = C' F (d t - d 2 )h + & E d 2 h 



where g is given by Eq. ( 96 1 . Within this approximation 



the metric tensor in Eq. ( 53 ) reduces to that of a planar 



geometry which allows us to identify the set of curvilinear 
coordinates s with the set x of Euclidean coordinates. 



In order to finally transform the equations ( 101 1 to the 



laboratory frame and to physical units, we define the 
interface position 



Z(x,f) = Vt + h(-K,t) , 



(102) 
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and the soft-mode concentration field 

C(r,i) ee C(r,t)-d[x,z- Z(x,t),t] (103) 
= C F [z- Z(x, t)} + c(x, z - Vt, t) + 0{h 2 ) . 



Even without recourse to the phase-field description it is 
clear from the equilibrium condition SH/SC = that 



U(z-Z) = C E {z-Z) 



(107) 



Assembling the results rt74h, and (101), written in 



physical units, we obtain, remembering the definitions 



(57), and (58), the set of equations 



is identical to the equilibrium concentration profile of the 
solute component. The equations of motion for the field 
variables Z(x, t) and C(r, t) are obtained in the form 



dtZ = K{ F + ad 2 Z 



d t Z = A [F 



-kJ + dz(c-U[*B(.»-Z)])d t U[$ E (z-Z)]\ , 
dtC = DV 2 {C -U[$ E {z- Z)]} (104) 



SH\ 
Jz I 



d t C = V-D(z-Z)V 



1 SH 

k ~5C ' 



(108) 



which determines the uniform propagation of the planar 
solidification front, including its stability behavior. 



Eqs. ( 104 ) present the requested interface description 



which can again be justified without going back to the 
phase-field description. In fact, in thermal equilibrium, 
where, instead of the driving force F, conveniently chosen 
Langevin forces enter [2B], the form of the equations ( 108 ) 



of the solidification process which generally holds at low 
growth velocities, but also applies to the rapid-growth 
regime of dilute metallic alloys. It, moreover, reconciles 
with the presence of a finite interface width, entering via 
the static phase-field profile <E> F . Remarkably, the latter 
is the only remnant of the initial phase-field description. 
This feature is a consequence of our expansion scheme 
which has been established by studying the phase-field 



essentially follows from the principle of detailed balance. 
As another aspect, we mention that, after insertion of the 



Hamiltonian (106), the equations of motion (108) are of 
a hydrodynamic type where capillary waves are, in the 
simplest-possible way, coupled to a bulk-diffusion field. 

Within this approach the uniform motion of a planar 
solidification front is, in dimensionless form, described by 
the unchanged first equation in Eqs. ( 74 1 . However, the 



model with the specific choices ph and pi. Once, this dra § force is now determined by the expressions 



scheme has been accepted, the interface model essentially 
follows by projecting the phase-field equation of motion 
to the soft-mode component $>' E . Application of the same 
procedure to the more general phase-field model, given 
by Eqs. ([!]), ( fl4| ), and (16), leads to an interface version 
where the quantities &e(z — Z) and D are replaced in 
proper positions by the functions 



G F (v) = - 1 l dCU'iOCpiO 

c F (C) -- 



(109) 



U(Z-Z) EE 

D(z - Z) = 



U[$ E {z-Z)] 
V[$ E {z-Z)\ 



where D(z — Z) has been re-scaled according to the last 
of the Eqs. The expression for C F (C) in Eqs. ( p9| ) 

derives from the once-integrated equation of motion 



(105) 



Similar input functions enter the description of grain- 
boundary motion by Cahn |19j who, however, uses an 
ideal-gas picture of the impurities in a Fokker-Planck 
representation of the system. In the following we are 
going to discuss some properties and applications of our 
generalized interface model. 



CAPILLARY- WAVE DESCRIPTION 



(110) 



which for the drag force G F (v) implies the representation 



(V ;( r) = - / dC^lC^t)} 2 . (Ill) 



From Eqs. (109) and (111) we conclude that generally 



Our first observation is that the new interface model 
can be written in terms of an effective Hamiltonian 



U = - / d 2 x{dZf 



d 3 r 



C-U(z-Z) 



(106) 



G F (0) = -2 7 , G F (oo) = , G F (v)<0 (112) 

where in the first equality we have used the behavior 
U(+oo) — U(—oo) = 2, following after re-scaling from 
Eqs. (|3|. Explicit evaluations of the quantities (109), of 
course, require specific choices of the functions U(z — Z) 
and D(z - Z). 
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For models, characterized by an arbitrary equilibrium- 
concentration profile U(z — Z), but a uniform diffusion 
coefficient D(z — Z) = 1, it is possible to determine all 
unstable eigenmodes of the solidification front. In order 



to demonstrate this, we start from the equations (96) 



and (101 1 with $ F (C) replaced by U'((). If, moreover, 



we consider perturbations of the form 

h(jc, t) = h(q, ui) exp (iq • x + uit) , (H3) 
C (C, x , t) = c(C, q, w) exp (iq ■ x + wt) , 

we encounter the set of equations 



(w +pg 2 )Mq,w) = -p7 / rfuJ7'(C)c(C,q,w) , 

*/ — OO 

( W + (7 2 -^ c -9 c 2 )c(C,q,^) = 
G F (C) (w + 9 2 )^(q, «) - C/'(C) <? 2 ^(q, w) 



(114) 

The last equation is a differential equation for c which, 
by the substitution 

p(C,q,^^^-C F (C) (115) 
can be converted into the more convenient form 

(co + q 2 -t-9 c -a 2 )p(C,q,c) = 
C/"'(C)-<Z 2 C/'(C) (116) 

where the source term on the right-hand side is directly 



expressed in terms of the input function U^Q). Eq. (116) 
has the solution 



p(C,q,w) - -U'(0 

X 2 -q 2 r°° 
A — n 



(117) 

dC'C/'(C')exp[A(C-C')] 



C 



+ T <t/'(C')ex P [MC-C')] 



X — fi 

with the characteristic roots A, given by 



A = 



i2 = 



-fl — V 



(118) 



Insertion of this solution into the first line of Eqs. (114) 
leads to the eigenvalue equation 

[u+pq 2 +E(q,w)]ft(cbw) = (119) 
where we have introduced a kind of self energy, 



/+oo 
1 
-OO 



duu'(0 [p(C,q,^) + c F (0] 



- P v[G F {v + X) -G F {v)] 



-P 



2X 



[G F (?; + A)+G F (A)] . (120) 



The final expression for S(q, w) has been obtained from 
Eqs. (62) and (117), and is identical to that found in 
Ref. 



|36j for the impurity-controlled motion of general 
domain boundaries. 



According to Eq. (120), non-trivial solutions h(q, uj) 
of Eq. (119) only exist under the condition 



+ pq 2 - pv[G F (v + A) - G F (v)} 



-P 



2A 



[G F (v + A) + G F (A)] =0 . (121) 



This relation determines the amplification rates ui(q) of 
all unstable eigenmodes of the model. It has a universal 
character, since it applies to a whole family of models 
with a globally uniform diffusion constant D, but with 
different equilibrium concentration profiles £/(£)• Given 
a specific form of such a profile, the explicit evaluation 
of Eq. ( 121 ) only requires the knowledge of G F (v) which 



follows from an analysis of the related one-dimensional 
growth scenario. 



By an expansion of the dispersion relation (121) in q 



and co we find two branches 0J\(q) and 0J2(q) with the 
behavior u>x(Q) = wi(0) = u}' 2 (0) = 0, and 



W2(0) 



G' F (v)-[G F (v)+G F (0)}/v-l 
l/p-G' F (v) 
±/p-G' f (v) 
{[G F (v) + G F (0)]M" ' 



(122) 



These expressions have a universal form in the same sense 
as Eq. (121). As demonstrated in Ref. [37], they even 
apply to models where the diffusion coefficient D(z — Z) 
has a constant value D in the liquid and the interface 
region, but is zero in the solid phase. Since for v — > 
the numerator in w"(0) shows the behavior 2|G F (0)|/i>, 
but approaches the value — 1 for v — > oo, the mode 0Ji(q) 
is unstable at low and stable at large velocities, provided 
the denominator in u>'{(0) is positive. Due to the first line 



in Eq. (74) this denominator is given by F'(v) which, as 
demonstrated below Eq. ( [77] ), can become negative. In 
this case the second mode 102(1) becomes unstable, since 
the denominator in W2(0) turns out to be negative in all 
applications of interest. An unstable point of the type 
q = 0, 0^2(0) > has previously been discovered by Cahn 
in a one-dimensional model of impurity-controlled grain- 
boundary motion [19] . in view of which we will denote 
the instability of the mode ^2(9) as a Cahn instability. 
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FIG. 8: Dispersion curve of the unstable Mullins-Sekerka-like 
mode oii(g), and of the stable mode 02(g) for the parameter 
values p = 100, 7 = 0.001, and v = 0.0036. 

10 5 co 




FIG. 9: Dispersion curve of the unstable Cahn mode 012(g), 
and of the stable mode wi (g) for the parameter values p = 100, 
7 = 0.01, v = 0.03. 




10 q 



FIG. 10: Dispersion curves of both unstable modes o>i(g) and 
012(g), arising for the parameter values p — 100, 7 = 0.01, 
v = 0.05. 



superposition of both instabilities. The complete loop in 
Fig. 10 shrinks to the point q = uj = 0, when the stability 
limit F'(v) = is approached, and it opens again in 
the stable region u) < 0, when F'(v) grows from zero 
to positive values. Qualitatively, the same behavior has 
been found by Braun et al. [38] in a phase-field model of 
solidification. 

A convenience of the capillary- wave description ( |106 ) 
and (108) is that it intuitively invites for applications 



to toy models with piece-wise linear functions J7(C) and 
-D(C)- A similar strategy has been used in the context 
of grain-boundary motion by Cahn |19) . and later, more 
extensively, by Hillert [33J . In the present case of binary- 
alloy solidification all consecutively considered families of 
such models have the equilibrium-concentration profile 



If, in a first application, we reconsider the model with 



p(0 = MC) + i . 0(0 = 1 



(123) 



we find, in accordance with the first line in Eqs. (76), 



G F (v) 



v + 2 

(v + iy 



(124) 



u(0 = e(c + 5)&(s-o 



C + S 



29(C-<5). (125) 



As shown in Fig. 11, this expression linearly interpolates 
between the two bulk phases, and for 5 = 1 mimics the 
the preceding model U(() — + l whereas for 6 — > 

it approaches the sharp-interface model U(() — 29(C)- If 
Eq. (125) is attended by a uniform diffusion coefficient 
D(Q = 1, Eqs. (109) lead to the expression 



G F (v) = 7[1 - 2v6 - exp (-2v8)] 



1 



(vsy 



and, as a consequence, to the the behavior 



(126) 



With that, a numerical evaluation of Eq. ( 121 ) in the 



unstable regime of both modes leads to the dispersion 
curves, shown in Figs. 8, 9, and 10 for real-valued rates 
loi, and u>2- They consecutively refer to the Mullins- 
Sekerka-like instability, the Cahn instability, and to the 



no) 



(127) 



11 3 

P Pc ' Pc ~ 2 £7 

From this result we see that the Cahn anomaly F'(0) < 
exists for finite values of 5, but disappears in the sharp- 
interface approximation. This is obviously the reason 
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FIG. 11: The potential U(() which, with a free parameter S, 
linearly interpolates between the liquid and solid phases. The 
limit 8 — > defines the sharp-interface limit in our approach. 



has the, again linearly interpolating, solution 



c F (0 = e(c + s)e(s-o ^ c + s (isi) 

1 + 2vo 6 



l + 2v8 



exp [-u(C - 6)} 



In view of the special value Cp (0) = 1/(1 + 2v5) we 
observe that the partition coefficient K(v) again has the 



form (80), however, with the new reference velocity 



V* = 



1_ Vd 
26 K K 



(132) 



In the measurements |34j of K(v) the quantity 8 may 
be used as a fitting parameter. Moreover, from the first 



formula in Eqs. (109) and the result (131) we obtain 



why the anomaly is even not mentioned in discussions, 
based on sharp-interface models, supported by local- 
equilibrium boundary conditions for the concentration 
field. As mentioned below Eq. (77), the anomaly is not 



excluded in the Aziz-Boettinger approach [2] where an 
internal structure of the interface region was effectively 
taken into account. 



In Ref. [37] we have considered a set of model systems 



which included the profile (125), but was complemented 
by a diffusion coefficient 



D(c) = e(( + 6) 



(128) 



Whereas for finite S the Cahn anomaly again appeared, 
we recovered in the limit 5 — > the instability, discussed 
by Misbah et al. [39] which is a kind of extension of the 
Mullins-Sekerka instability [3] from the diffusion- into the 
kinetics-limited regime. 

As a final application we now consider a class of models 



where the equilibrium concentration (125) is attended by 
the diffusion coefficient 



D(C) = o U(() 



(129) 



which, consequently, also linearly interpolates between 
the two bulk phases. This behavior mimics the result 
of numerical calculations, derived in Refs. |40| . and [41] 
from the mean-square particle displacements. Adopting 



the form (129), one easily verifies that the steady-state 
equation 



c'AO + n^friO = u '(0 (iso) 



G F (v) 



l + 2v8 



(133) 



which, via the force-balance relation in Eqs. ( 74 ) , implies 



F(v) = 1--' 



1 



p p c 1 + 2vS 



Pc 



1 

4^7 



(134) 



This function is convex in the sense F"(v) > 0, and has 
a minimum at 



(135) 




F„ 



so that the Cahn anomaly F'(v) < occurs in the whole 
regime v < v m . By means of Eq. (98) the value F m 
corresponds to the temperature 




(T L - T s ) 



(136) 



where all temperatures refer to the density Cs- The value 
T m limits the temperature range, up to which a meta- 
stable kinetics-limited solidification is possible inside the 
diffusion-limited regime. In the constant-miscibility-gap 
approximation the line T m (C) is, due to Eq. (fsb , located 
below the Baker-Cahn line To(C), and approaches it in 
the limit p — > oo. 



If Eq. ( 134 ) is rewritten in terms of the scaled variables 



1 V 

w — 

it assumes the simple form 



. / 



27 



-F 



(137) 
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arises in case of the anomaly F'(v) < 0. A regime with 
this behavior occurs in all considered models, except in 
the sharp-interface limit where only the Mullins-Sekerka 
instability survives. The behavior F'(v) < even exists 
at the origin v — which has also been noticed in many 
approaches, based on a phase-field description. 

In view of the anomalous behavior F'(0) < one may 
question the consistency of our model, and also that of 
standard phase-field models, with the basic principle of a 
positive entropy production. It has been shown, however, 
by Bi and Sekerka that, for a general class of phase- 
field models, the entropy production is positive. Applied 
to our model, and using our notations, their expression 
for the entropy production can be written in the form 



FIG. 12: Kinetic spinodal line w(f), enclosing the region 
where the Cahn instability f'(w) < occurs in the /, w-plane. 
The limiting point w(— 1) = corresponds to the adiabatic 
limit in the one-component system. Also shown are some 
curves f{w) for different values of the parameter p/p c - 



(138) 



where only the ratio p/p c enters as a tunable parameter. 
Elimination of this parameter from Eqs. ( 135 ) leads to 
the trajectory 



iv = V-f i- V-f 



(139) 



in Fig. 12, enclosing a regime of unstable behavior in the 
sense f'(w) < 0. The outside region, up to the line / = 0, 
is a a regime of meta-stable solidification. Accordingly, 
the trajectory may be considered as a spinodal line of 
kinetic origin which starts at the origin w = / = for 
p — p c , and approaches the point w — 0, / = —1 for 
p — > oo, corresponding to the limit Dl — > 0. 



Remarkably, the result (1391 precisely agrees with an 



analogous stability limit, derived by Umantsev for the 
non-isothermal solidification of a pure one-component 
substance 22]. Since, in this case, the diffusion constant 
Dl has to be replaced by the heat-diffusion coefficient 
Dt, the above-mentioned limit Dt — > corresponds to 
the process of adiabatic solidification. 



DISCUSSION 

The most remarkable implication of our approach is 
the existence of the Cahn instability ^2(9) > which 



n= / d 



V 



sc 



(140) 



which obviously is positive. In order to ensure that the 
expansion scheme, underlying our capillary-wave model, 



does not endanger this property, we apply Eq. ( 140 1 to 



the one-dimensional steady-state solidification, replacing 
everywhere $_f(C) by the static profile $>e(0- The result 
for the entropy production per unit area then reads in 
dimcnsionlcss units 



V 



{d^Ef+D^C'p-U'f 



G F (v) 



v[F(v) - G F (0)] (141) 



where we have used the expression (26) for the 
representation of Gf(v) in the first line of Eqs. (109), 
the relation (110), and the force balance in Eqs. (74). 



In view of Eq. (112) the quantity (141 ) again is positive 



which remains true in the expanded form 



7T « v[2j + vF'(0)} 



(142) 



and, consequently, demonstrates that the strange-looking 
behavior F'(0) < is not in conflict with the second law 
of thermodynamics. 

To complete our presentation, we add a few estimates 
which border the range of applicability of our approach. 
First of all, it is clear that the two small parameters 7 
and v/p of our approach can be tuned by the miscibility 
gap AC and by the driving force F, respectively. From 
Eqs. (10) and (20) we know that 



7 



£L ( dC L 
4ct ^ dT 



AC 

Tm 



(143) 



According to Turnbull 43 there is a correlation between 
the surface tension and the latent heat of the form 



a = Ct L a 



(144) 
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where a measures the average atomic distance, and, for 
essentially all metals, Ct ~ 0.45. Assuming £ w 1.8 a, 
and expressing the miscibility gap by the temperature 
gap 



AT=T L (C s )-T s (Cs) , 



Eq. ( 143 1 reduces to 



7 



AT 



(145) 



(146) 



For the values AT = 10 K and T M = 1000 K this yields 
7 R3 0.01 which is sufficiently small, and identical to the 
value, underlying Figs. 9 and 10. The value p = 100, also 
used in these figures, has been adopted from Ref. [2"T] . 
and, in the rapid growth regime v « 1 implies v/p ~ 0.01 
for our second expansion parameter. 

As a final remark we point out that the capillary-wave 



model, defined by Eqs. ( 106 1 and ( 108 1, is amenable to a 



variety of generalizations without going back to a phase- 
field description. One generalization in the surface part 



of the effective Hamiltonian ( 106 1 is the replacement of 



the integral expression by the exact area of the interface 
which in the equation of motion leads to the appearance 
of the full mean curvature. Simultaneously, one may also 
incorporate an anisotropy of the surface tension. Further 
generalizations are the inclusion of the energy density or 
other bulk fields which then generates a number of cross 



couplings in the Hamiltonian. Parallel to this, Eqs. ( 108 ) 



will be replaced by an enlarged set of equations of motion 
which in general is equipped with a matrix of Onsager 
coefficients. A consistent discussion of the directional 
solidification of a dilute binary alloy requires a coupling 
of the energy density to the interface position which we 
will consider in a forthcoming paper. 

We are grateful to Efim Brener for calling our attention 
to the consistency problem, examined in our discussion. 
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of Diisseldorf for its warm hospitality. This work has 
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